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1. INTRODUCTION 


The source of the broad-scale, slowly varying, large-amplitude 
portion of the geomagnetic field observed at and near Earth's surface is 
widely held to be electric current flowing in Earth's electrically 
conducting core and, to a lesser extent, electrically resistive mantle 
and magnetically permeable crust. Maintenance of this current against 
resistive dissipation is usually attributed to motional induction in the 
liquid outer core: motion of the magnetized fluid relative to the solid 
mantle and crust generates electromotive force and thus electric current 
(LARMOR, 1919). If the curl of this electromotive force field is 
neither zero nor perfectly balanced by magnetic diffusion, then the 
fluid motion will induce slow or secular change of the geomagnetic field 
in accord with Faraday's law (ELSASSER, 1946a, b, 1947) . Indeed, if an 
initial magnetic field, the relevant material properties, and the 
evolving fluid velocity field were all known within the Earth, then a 
straightforward, classically deterministic prediction of geomagnetic 
secular change could be made, The inverse problem of estimating the 
fluid motion within the Earth given real observations of the magnetic 
field near its surface and rough estimates of key material properties 
is, of course, not quite so straightforward. 

Both the rich variety of geomagnet i ca 1 ly significant Earth 
properties and processes and the mathematical subtlety of the implied 
geophysics present formidable obstacles to core flow estimation. It has 
thus seemed appropriate to seek a qualitative understanding of 
geomagnetic secular variation in terms of a simple model, attempt 
quantitative explication of the observations, and to then systematically 
restore essential geophysical detail. In addition to using a simple 
magnetic earth model, I have further restricted my study to the 
geomagnetic effects of induction by hypothetically steady fluid motion 
(and steady magnetic diffusion) near the top of Earth's core. 

The purpose of this series of papers is to share the findings of my 
decade-long investigation of steady induction effects in geomagnetism. 
These findings include qualitative geophysical arguments, rigorous 
mathematical proofs, technical development and application of new 
methods, and quantitative numerical, results . Considered space is 
devoted to the technical introduction, background theory, mathematics, 
and physics needed to make these findings clear to the general reader. 
Indeed, some reviewers have demanded more detail precisely where others 
have demanded less. I hope that experts will share my belief that the 
findings transcend occasionally divergent stylistic preferences. 

1.1 A Simple Magnetic Earth Model and Some Enabling Conditions 

Following ROBERTS a SCOTT (1965), BACKUS (1968), BENTON (1979) and 
others (see, e.g., HIDE & MALIN, 1981; VOORHIES, 1984, 1986a, 1987a; 
BLOXHAM & JACKSON, 1991), attention is focused upon the fluid motion 
near the top of the core by adopting a simple model of the magnetic 
Earth: the source-free mant le/f rozen-f lux core model (or SFM/FFC model) . 
In this model, a rigid, impenetrable, electrically insulating mantle of 
uniform magnetic permeability surrounds a spherical, inviscid, perfectly 
conducting outer core in anelastic flow. The collective use of these 
eight suppositions to account for recent broad-scale geomagnetic secular 
change is expected to yield errors of about 7%-as will be discussed in 
Part II (VOORHIES, 1986c, 1987c) . In the SFM/FFC model, geomagnetic 
lines of force that thread the core-mantle boundary (CMB) are rooted in 
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the fluid at the top of the core; lateral motion of this fluid induces 
secular change by advecting the two footpoints of each such magnetic 
field line so as to reconfigure the scaloidal magnetic field outside the 
core. In and just above a SFM, the magnetic flux density vector B is 
the negative gradient of the scalar potential V: B = -Vv. 

Granting these simplifying suppositions, even complete and perfect 
information about the geomagnetic field at Earth's surface does not 
allow unambiguous determination of the fluid motion at the top of the 
core (ROBERTS & SCOTT, 1965; BACKUS, 1968). The radial component of the 
induction equation at the top of a frozen-flux core attributes the time 
rate of change of the radial magnetic flux density component to the 
surface convergence of the product of this component with the lateral 
fluid velocity ( 3 1 B r = _ ^ s * B r v s )# so an Y cryptic flow which induces no 
secular change is not detected (any v ' such that V g « B r v s ' = 0). 
Despite this toroidal ambiguity in B r v s (BACKUS, 1982), some components 
of this motion could be so determined at special locations {BENTON, 
1981). Moreover, this motion could be uniquely determined from such 
information under certain conditions: e.g., if and when the motion is 
steady (VOORHIES & BACKUS, 1985) or, in broad areas, if and when it is 
tangentially geostrophic (HILLS, 1979; BACKUS & LEMOUEL, 1986). 

Such enabling conditions, and the SFM/FFC model itself, can be 
viewed as approximations to a more realistic core geodynamo model: 
supposition of steady flow is accurate to lowest order in a temporal 
Taylor series expansion of the fluid motion; tangential geostrophy might 
hold to lowest order in core dynamics.' The inaccuracy of these 
approximations might be demonstrated quantitatively by analysis of 
geomagnetic observations. If only to establish the importance of other 
effects, such conditions may then be treated as hypotheses to be tested 
against geomagnetic data (or models thereof) in the context of the 
SFM/FFC model (or refinements thereof) . 

Geomagnetic data will be neither complete nor perfect within the 
forseeable future, so geomagnet ic~ test ing of such hypotheses is further 
embedded in at least one additional supposition regarding the 
geomagnetic field itself. For example, an initial condition is needed 
to test the steady motions hypothesis; therefore, the rigor of 
geomagnetic tests of this hypothesis is limited by the completeness and 
accuracy with which an initial geomagnetic field can be specified. 

Magnetic field measurements do not isolate the slowly varying, 
broad-scale portion of the field-the "core field" of interest here. 
This isolation is imperfectly achieved with truncated spherical harmonic 
models derived to fit select, often low-pass space- and time- filtered, 
geomagnetic data (see, e.g., LANGEL, 1987); however, such models do 
represent the relevant portion of the observations. In recognition of 
the effects of external and non-core internal sources, the use of such 
models is here preferred to the use of raw data. Spherical harmonic 
models are convenient to work with; yet substitution of spatially 
complete, but spectrally incomplete, truncated spherical harmonic models 
for (in effect) spectrally complete, but spatially incomplete, 
measurements is not without pitfalls. For example, spherical harmonic 
(Gauss) coefficients of the scalar geomagnetic potential are uncertain 
and cross-correlated (LANGEL, ESTES & SABAKA, 1989; LANGEL, 1991); 
moreover, coefficients that are not estimated are not necessarily 
assumed to be zero. Spherical harmonic models may represent reduced, 
filtered, and analyzed data; but neither derived model parameters such 
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as Gauss coefficients nor annual averages of magnetic observatory data 
should be confused with the result of a physical measurement. 

1.2 Some Motiv ation for Considering Steady Core Surface Flow ^ ^ 

The supposition of statistically steady surficial core flow during 
intervals of a few decades or more is adopted for a variety of reasons. 
Theoretical arguments (invoked a posteriori) based upon ties between core 
surface flow and extremely slowly varying thermal, topographic, and 
compositional anomalies in the deep mantle support this supposition. 
For example, obstacles to a persistent westward flow posed by 
depressions of the CMB could result in an effectively steady, yet 
spatially complicated flow near the core surface. RUFF & ANDERSON 
(1980) and BLOXHAM & GUBBINS (1987) suggest "thermal core-mantle 
interactions," whereby laterally varying heat flow from the core to the 
thermally heterogeneous deep mantle establishes a steady, superficial 
core circulation. JAULT & LEMOUEL (1989) suggest "topographic locking," 
whereby the pattern of perturbation pressure driving core surface 
circulation becomes locked to the topography of the CMB. Another 
possibility is "Lorentz linkage, * whereby strong Lorentz forces inhibit 
the motion of fluid parcels lying on field lines which thread regions of 
anomalously high conductivity in the deep mantle. Such patches of 
relatively stagnant fluid may be separated by mobile streams wherein 
fluid parcels are but loosely coupled to overlying areas of weak mantle 
conductivity. Such streams could slowly but steadily transport fluid 
laterally from large-scale source regions towards sinks. The sources 
could be buoyant FeO or FeS enriched reservoirs trapped near the CMB 
after columns or plumes from the inner core boundary have impinged upon 
the CMB (suggested by S. I. BRAGINSKY, pers. comm., 1989); the sinks 
might be compositionally depleted due to wOstite underplating of the 
initially low-conductivity region overlying the stream. One might also 
imagine steady drift of an otherwise standing wave field. 

Facts, practical considerations, and empirical arguments may offer 
more compelling reasons. For example, supposition of surficially steady 
flow is simple, annihilates a key formal ambiguity (VOORHIES & BACKUS, 
1985) , and retains the lowest order term in a temporal Taylor (or 
Fourier) series expansion of the fJLow which, by definition, always 
merits consideration. BLOXHAM (1987a) suggested that all geomagnetic 
estimates of core surface flow suppose steadiness on some time scale. 
Those which do not (e.g., BENTON & CELAYA, 1991) still suppose some 

other form for the time dependence and make at least one quasi-steady 

approximation. Moreover, the primarily kinematical supposition of 
steady flow provides an alternative to specific dynamical assumptions 
and might thus shed some light on core dynamics-a view strongly 
advocated by WHALER (1991). It also brings all relevant data to bear on 
the problem of estimating a single core surface flow-a statistically 
interesting property in light of practical ambiguities resulting from 
the finitude of geomagnetic measurements, each of which is influenced by 
myriad geomagnetic phenomena. Perhaps most importantly, it turns out 

that such a single steady core surface flow explicates recent secular 

change to better than first-order accuracy. 

Most of the recent secular variation (SV) indicated by broad-scale 
spherical harmonic models of the observed geomagnetic field can be 
explained quantitatively in the SFM/FFC model by a steady surficial core 
flow (VOORHIES & BENTON, 1983; VOORHIES, 1984, 1986a)-especially when 
the derived flow is not artificially constrained to induce minimal 
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narrow-scale SV (VOORHIES, 1986b). Published studies by K. Whaler, J. 
Bloxham, and others, and unpublished, if not wholly reported, studies by 
this author and others confirm this fact. In particular, this fact 
remains when the steady flow is constrained to be, in effect, 
surficially geostrophic (VOORHIES, 1986d; 1991); when the spatial 
complexity of the flow is damped (BLOXHAM, 1986, 1987a, b, 1988a, 1989; 
VOORHIES, 1987b, 1988a; WHALER & CLARKE, 1988); and when both 
geostrophic and smoothness constraints are enforced (VOORHIES, 1987c, 
1988a, 1989, 1 99 1 ; BLOXHAM, 1987c, 1988b). This fact supports the 
argument that time-varying flow near the core surface and departures 
from the SFM/FFC model in the real Earth are of secondary importance in 
the correct explication of recent SV. It also supports the argument 
that the flow near the top of the core is similar to the estimated flows 
-many estimates derived by different groups using diverse methods and 
models being positively correlated (WHALER, 1990) . 

These arguments seem provocative enough to motivate an in-depth 
exposition of both the physical foundation upon which, and the methods 
by which, this fact is established. Moreover, colleagues suggest that 
the approach, theory, method, techniques, and results of some recent 
attempts at steady surficial core flow estimation merit more than 
mention at meetings. There are, of course, other reasons to investigate 
steady induction effects in geomagnetism-one being to learn to 
distinguish such effects from other effects. 

Limited quantitative success of a simple interpretation does not 
prove it correct. Indeed, interpretation of recent SV solely in terms 
of the SFM/FFC model and steady surficial core flow is in error due to 
the extreme simplicity of these seemingly relevant idealizations. 
Nonetheless, the demonstration that these idealizations provide a 
powerful tool for quantitative explication of SV reflects favorably upon 
the power of any more realistic core geodynamo model. It also shows 
that statistically steady surficial core flow can dominate geomagnetic 
field behavior over appreciable intervals, may thus be a component of a 
more realistic core geodynamo model, and might even be a minor feature 
of an accurate core geodynamo theory. The geomagnetic effects of steady 
motional induction thus merit study not only to establish interesting 
facts, estimate core surface flow, and provide a fascinating contrast to 
other possible sources of SV, but because an understanding of these 
effects is important, and perhaps essential, to the correct 
interpretation of SV— if not the phenomenon of geomagnetism itself. 

1.3 Outline of the Series on Steady Induction Effects 

To investigate steady induction effects in geomagnetism, and to 
help establish the importance of other effects, it is often useful to 
regard the SFM/FFC model as a first approximation and to treat the 
supposition of steady surficial core flow as a hypothesis. To test 
hypotheses against observations it seems appropriate to (a) understand 
both; (b) develop a satisfactory method for mimicking the relevant data 
in accord with the hypotheses; (c) apply the method to make quantitative 
predictions; and (d) subtract predicted from observed values and measure 
such residuals in units of the estimated uncertainty in the data. This 
series of papers describes an effort to meet these requirements using 
relevant geomagnetic field models to test the SFM/FFC earth model (and 
refinements thereof) and the hypothesis of (piecewise, statistically) 
steady surficial core flow. 
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In Part I , the steady surficial core flow estimation problem is 
examined and solved in the context of the simple SFM/FFC earth model. 
The present, introductory paper (IA) reviews the basic theory and 
develops some implications of the steady motions hypothesis. Paper IB 
develops a method for solving the non-linear inverse motional induction 
problem posed by the hypothesis of (piecewise, statistically) steady 
core surface flow and the adoption of an initial geomagnetic condition. 
This inverse problem is generally non-linear because of the nature of 
the solution to the forward problem and because neither the models nor 
the observations they represent are either complete or perfect. Paper 
IC describes application of this method to the Definitive Geomagnetic 
Reference Field (DGRF) models (IAGA, 1988). Paper ID presents numerical 
results of applying the method and conclusions drawn therefrom. 

In Part II, the SFM/FFC model is reexamined; errors induced by 
oversimplifying suppositions are assessed and targeted for systematic 
elimination. In Part III, the supposition of perfect core conductivity 
is replaced with that of steady magnetic flux diffusion near the top of 
a resistive core and an effort to allow for, and indeed estimate, deep 
mantle electrical conductivity is described. 


2. THEORY 


Many geomagnetic estimates of core surface motions rely on the 
magnetic induction equation for a fluid medium of isotropic magnetic 
permeability |l and steady, isotropic electrical conductivity G 

3 t B = Vx(vxB) + Vxo" 1 (Vx(|i _1 B) ] < la ) 

where B represents the magnetic flux density vector, dj-B its partial 
derivative with respect to time t, and v the fluid velocity vector. The 
first term on the right of (la) describes motional induction; the second 
describes magnetic flux diffusion. The optional appendices offer a 
derivation of this equation for general readers and an analysis of 
conditions under which it is a useful summary of the Maxwell equations. 
For homogeneous [i and uniform O (la) reduces to 


d t B = Vx(vxB) 


T|V 2 B 


(lb) 


where T| s [|IC]”^ represents the magnetic diffusivity and the vector 
identity VxVxA = VV*A - V 2 A has been applied to solenoidal B. 

In spherical polar coordinates (r,0,<|>) with orthonormal unit 
vectors (r,0,$), let A represent a vector field with radial component A r 
is r*A and surficial component A s = A - r A r . Let V g » denote the surface 
divergence operator on the sphere of ^radius r: V S *A = V*A - r 9 r r A r = 
[rsinO] -1 [SgAgsinG + . Because r* [Vx(AxZ) ] = [A r Z g - Z r A s ], the 
radial component of (lb) is 


a t B r + V s *[B r v s - UB S ] = Tlr _ 1 V 2 rB r = d r (1c) 

where B r and u = v r represent, respectively, the radial components of B 
and v, B s and v s the surficial portions of B and v, and d denotes the 
radial maqnetic flux density diffusion. Equation (lc) is the radial 
induction equation for uniform, isotropic media. To restore effects of 
displacement currents, add Xg3 t 2 B r to the left of (lc) - where Xg = e /o 
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is the charge relaxation time and £ is the uniform, isotropic dielectric 
permittivity (Appendix A); omission of this term filters out 
electromagnetic radiation in accord with the quasi-steady approximation. 

2.1 Radial Induction at the Top of a Spherical Core 

The magnetic flux density vector B(r,t) and the fluid velocity 
vector v(r,t) within the Earth may vary with time t and position vector 
r reckoned in geobarycentric spherical polar coordinates (r,9,<|>) fixed 
to the solid Earth; henceforth r is radius, 0 is colatitude south from 
the reference pole of rotation, and <)> is longitude east from the 
Greenwich meridian. The upper part of the outer core is treated as a 
low-viscosity liquid of uniform |i and uniform a wherein (1) applies. 
The CMB is approximated by the sphere of radius b = 3480 km-the mean 
PREM core radius (DZIEWONSKI & ANDERSON, 1981) . The surrounding mantle 
is treated as a perfectly rigid, impenetrable, electrical insulator of 
uniform magnetic permeability wherein, perforce, v = 0 and scaloidal B = 
-W. The viscous boundary layer just beneath the CMB is treated simply 
as a sheet vortex; points on the spherical surface just beneath it lie 
at the top of the free-streaming core and are assigned positions r = b. 
The kinematic boundary condition is applied at b; thus v(b,t)*r s u(b,t) 
= 0, v(b,t) = v g (b,t) is tangent to b, and (lc) at (b,t) is 

3 t B r (b,t) + V s *[B r (b, t) v s (b, t) ] = T) [r _1 V 2 rB r * (r, t ) ] s d r (b,t) (2) 

The asterisk on B r * emphasizes that only portions of the field that are 
not purely scaloidal contribute to d r . 

Across a spherical boundary B r is continuous, 3 r B r can jump with B s 
due to a sheet electric current, and 0 r 2 B r will jump with the lateral 
electric current density J s (see Appendix A) . The radial component of 
any part of the field which is not purely scaloidal, say B r , vanishes 
at the spherical base of a source-free mantle, so it is negligible at b 
by continuity; however, d r 2 B r * will typically be non-trivial at b, as 
will d r (b,t) . When the core is treated as a perfect conductor T| 
vanishes and d r (b,t) is zero-unless a sheet current and jump in d r B r 
dictate otherwise. Though radial flux diffusion at the top of the free- 
stream may well account for about 6% of broad-scale secular change 
(VOORHIES, 1986c; 1988 unpublished, 1989), d r (b,t) is omitted in part I. 
The result is the celebrated ROBERTS & SCOTT (1965) equation 

5 t B r (b,t) + V s *[B r (b, t) v s (b, t) ] = 0 (3) 

which describes frozen-flux motional induction at the surface of a free- 
streaming spherical outer core. 

2.2 Kinematics 

Mass density p(r,t) is treated as the sum of a mean state p Q (r) a 
fluctuating perturbation p'(r,t). Mass conservation at (r,t) implies 

d t p' + V» [ p ' v ] + V*[p Q v] = 0 (4) 

where V* is the full divergence operator. The anelastic approximation 

V»[p 0 (r)v(r, t) ] = 0 (5) 

filters out high-frequency acoustic radiation in much the same way that 
omission of displacement currents from the Ampere -Maxwell law filters 
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out high-frequency electromagnetic radiation. It does not require 
vertical advection of the mean stratification and is considered more 
appropriate to the slow outer core flow of interest here than is the 
supposition of solenoidal flow (VOORHIES, 1987c, 1988 unpublished ms.) . 

Equation (5) differs from equation (1.36) of GUBBINS Sc ROBERTS (1987) in 
that advection and flow convergence are allowed to change p'(r,t). 

In a spherical outer core, p Q (r) = p Q (r) , so (5) and the kinematic 
boundary condition at b imply incompressible flow at b: V^vl^ = 0. Then 
V«v_(b,t) = -3 r ,u(b,t), surficial convergence (or confluence V *v < 0) 
implies downwelling (o r u > 0), and (3) reduces to 

d t B r (b,t) + v s (b,t)*V s B r (b,t) = B r (b,t)d r u(b,t) . (6) 

This is the usual frozen-flux radial induction equation at b. 

In this frozen-flux core ( FFC ) approximation, the mean square 
radial magnetic flux density linking the core changes only when fluid 
downwelling correlates with the squared radial magnetic flux density. 
To see this, let <\|/( r,t)> represent the mean value of any scalar field 
\|f(r,t) averaged over the sphere of radius r 

271 K 

<\y(r,t)> s ( 47t) _1 J J [\|/(r, t ) ] sin0d0d<j) (7) 

0 0 

and let \|/(r,t) rms denote the root mean square value of \|/{r,t) averaged 
over the sphere of radius r, <\|A(r , t ) 2>l/2 . Note that <V g ®A> = 0 for any 
vector field A(r,t) with single-valued differentiable components. Now 
consider the mean square radial field averaged over b 

27C k 

(47I) _1 J J [B r (b, t) ] 2 sin0d0d<]) = < [B r (b, t ) ] 2 > s [B r (b, t ) rms ] 2 . (8) 

0 0 

The value of (8) indicated by a finite set of Gauss coefficients is 
easily calculated in the SFM model (LOWES, 1966; VOORHIES, 1984) . With 
(3) and V s *v s (b,t) = -3 r u(b,t) the time derivative of (8) is 

<2B r (b,t)d t B r (b,t)> = <2B r (b,t) [B r (b f t)3 r u(b,t) - v s (b, t ) •V g B r (b, t ) ] > 

= <2B r 2 3 r u - V s .[B r 2 v s ] + B r 2 V s .v s > 

= <[B r (b,t) ] 2 a r u(b,t)>. (9) 

Indeed, at (b,t) we have 2<B r 3 t B r > = <B r 2 8 r u> = 2<B r v g * V g B r > ; yet the 
contribution to (9) from lateral advection is one half, but opposite, 
that from downwelling. 

In the special case of purely toroidal flow at the surface of a FFC 
v s (b,t) = v T ( b , t ) = [V g T(b, t) ]x [r] , where -T is the streamfunct ion ; then 
V s * v t ( b , t ) = 0, there is no downwelling, and the right-hand side of (6) 
and (9) vanish with 3 u(b,t). Any purely toroidal flow at the surface 
of a FFC therefore conserves the mean square radial flux density 
averaged over the surface. This is obviously so for rigid rotations, 
but also holds for arbitrarily complicated toroidal flows-al lowing a 
quick check on the compatibility of estimates of B r (b,t) with the 
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hypotheses of purely toroidal FFC surface flow. Alternately, if the 
flow is toroidal, then there is a second way to locate the core 
magnetically: seek the sphere upon which <B r (r,t) 2 > is unchanged over 
some time interval. (The first way is described by HIDE & MALIN (1981), 
VOORHIES Sc BENTON (1982), VOORHIES (1984), BENTON Sc VOORHIES (1987)). 
Mean square radial field conservation is also easy to impose on models 
of the evolving geomagnetic field; such biased models might be useful 
for seeking purely toroidal flows. Of course, a geomagnetic field model 
constrained to be smooth in a well-defined sense that also yields small 
values of (9) will, perhaps unintentionally, be more compatible with the 
no-upwelling hypothesis than models that are not so constrained. 

It is clear that supposition of sphericity is not needed (BACKUS, 
1968) to equate the time rate of change of the mean square magnetic flux 
density normal to the closed surface of a perfect fluid conductor 
bounded by a rigid, impenetrable exterior with the mean value of the 
product of the squared normal flux density and the surface convergence 
of the surficial fluid velocity averaged over the surface. 

2.3 Simplified Dynamics 

The dynamics of the outer core is a subject of debate; however, a 
reasonable approximation to the equation of fluid motion therein is 

p[3 t v + v*Vv + + 2Qxv] = V*FI + pg e + JxB (10) 

where Q is the bulk angular velocity of the solid earth; n is the 
hydrodynamic stress tensor (including pressure and bulk and shear 
viscous stresses) and V*n is its divergence; g e = - O.slflxrl 2 ] is 
the effective gravitational acceleration for Newtonian potential % and 
centrifugal potential -0.5|Qxr| 2 ; and J is the electric current density 
(from A3d) . Physically, V«n is the macroscopic representation of net 
microscopic electromagnetic force densities and the Lorentz force 
density JXB represents the divergence of macroscopic magnetic stresses 
(see, e.g., VOORHIES, 1991). 

The mean state of the outer core is taken to be the equilibrium 
between forces caused by hydrostatic pressure and by gravitational 
attraction and centrifugal effects: Vp Q (r) * p Q (r)g e (r), where V is the 
gradient operator and p Q is -1/3 the trace of diagonal n for a Newtonian 
fluid at rest in the reference frame rotating at This mean state is 
subtracted from (10). High-frequency fluctuations in Q and tidal 
effects are omitted (in Part I) . For a modest to weak magnetic field, 
scale analyses of the residual perturbation momentum equation in the 
upper part of the core suggest that, just beneath a thin (8 cm) viscous 
sub-layer, the relative, advective, and precessional pseudo-force 
densities on the left and viscous, centrifugal, and Lorentz effects on 
the right of (10) contribute little to the primarily geostrophic balance 
between Coriolis, perturbation pressure, and radial buoyancy forces 
(see, e.g., HILLS, 1979; LEMOUEL, 1984; BENTON, 1985; VOORHIES, 1991). 
Then with p' « p Q , omission of perturbations in g 0 , and g e = -rg, 


2p o n o xv + Vp ' + p f gr - 0 (11) 

where the prime indicates the (non-tidal) perturbation relative to the 
mean state. A modest magnetic field at the top of the free-stream is 
consistent with an effectively SFM, but need not imply a weak field at 
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depth (e.g., below^a 90-km-thick magnetic boundary layer). With £2(t) =* 
J2 o = Q o (rcos0 - §sin0), the radial component of the curl of (11) 
evaluated at b reduces to the geostrophic radial vorticity balance 

V 

3 r ucos0 + -sin0 » 0 (12) 

b 

which has long been used to constrain some numerical estimates of core 
surface flow. Physically, (12) requires downwelling to be accompanied 
by poleward flow except at the equator (where v = 0) and at the poles 
(where d r u = 0) (VOORHIES, 1987d) . Note that (12) approximates the full 
radial vorticity equation under somewhat less restrictive conditions 
than those under which (11) follows from (10) (BENTON, 1985; VOORHIES, 
1991) . Flows obeying (12) are "surficially geostrophic" and include 
tangentially geostrophic flows defined by BACKUS & LEMOUEL (1986) . 

3. STEADY MOTIONAL INDUCTION: THE FORWARD PROBLEM 

For steady flow, (3) reduces to 

a t B r (b,t) + V s *[B r (b, t)v s (b) ] = 0 (13a) 

For steady anelastic flow, (6) reduces to 

a t B r (b,t) + v s (b)*V s B r (b,t) = B r (b,t)3 r u(b) . (13b) 

Equations (13) also describe frozen-flux motional induction during an 
interval t Q <5 t < t^ when the flow is steady in the statistical sense: 
when correlations between fluctuations about short time averages of v s 
and B r contribute negligibly to short time averages of 3 t B r and the 
short time averages of v s do not vary during the long time interval 
[ t o , t £ ] (VOORHIES, 1986a, section 2.2). A flow that is statistically 
steady during consecutive intervals is piecewise statistically steady; 
though v s may change between intervals, (13) holds within each interval. 
Effectively piecewise steady flow might arise physically when long 
intervals of statistically steady flow, hence statistically balanced 
forces, are punctuated by rapid shifts to a new flow configuration. 

Given initial condition B r (b,t c ) and v s (b) the forward solution to 
(13) is of transcendental exponential operator form [Voorhies, 1986b] 


-< t - t o )V s #v s (b) 

B r (b,t) = (e }o{B(b,t 0 )} (14a) 

or, perhaps more clearly, 

-(t-t c )V s *[v(b) 

B r (b, t) = {e }oB(b,t 0 ) ] . (14b) 

This is a special case of the general solution to forward steady 
motional induction problems presented by HOYNG (1985) . To derive (14) 
time differentiate (13a) and back substitute k-1 times: 

a t 2 B r + v s *[v s a t B r ] = a t 2 B r - v s *[v s {V s »[v s B r i n = o (isa) 
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0 


(15b) 


d t 3e r + V s* tv s 3 t 2B r ] = 3 t 3e r + V f v s { V t v s < V ( v s B r ] 5 1 H = 
a t k B r = (-l) k {V s *[v s } k oB r ] (15c) 

where (15c) follows by logical induction. Equation (15c) at t Q is 
substituted into the Taylor series expansion for B r (b,t) about t Q 
oo 1 ^ 

B r (b,t) = z — (t-t Q ) a t k B r (b,t 0 ) 

k=0 k! 

oo (-l) k 

= 2) (t-t o ) k (V s .[v s ) k 0 B r (b,t o )] 

k=0 k! 

«> 1 

= Z -{-(t-t 0 )V s .[v s ) k oB r (b,t 0 )] 
k=0 k! 

~ ( t-t ) V • [ v (b) 

= {e )oB r (b / t 0 )] 

which is (14b) or, with the transcendental exponential operator taking 
precedence over the surface divergence operator, (14a). 

One feature of evolution equation (14) is that action of a given 
steady surface flow on different initial radial magnetic flux conditions 
induces predicted flux configurations which can diverge exponentially 
with time in some regions and converge exponentially with time in other 
regions. Consider the iterative map obtained by forward numerical 
solution of the steady motional induction problem (13b) with an 
arbitrarily small time step 5t 

B r (b,t+St) = B r (b,t) + O t B r (b, t) ]St. (16) 

With time measured in units of 5t 


B r (b,t+1) = B r (b,t) + d t B r (b,t) 

= B r (b,t) - v(b)»V s B r (b,t) 

dB r (b,t+l) 

= 1 - V s »v(b) = 1 + d r u 

dB r (b,t) 

and the Liapunov exponent is 

1 N dB r (b,t+l) 1 

L = lim - X In I — 1 = lim - 

N->ooN t = l dB r (b,t) N^ooN 

Restoration of the original measure 

L = In 1 1 + 5t3 r u| . 

The sign of the Liapunov exponent i: 


(17a) 

- B r (b f t)V s *v(b) (17b) 

(17c) 

N 

X lnll + 3 r u| = lnll + 3 r ul . (18) 

t = l 

of time yields 

( 19 ) 

the sign of the downwelling! 


! 
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For downwell ing 3 r u >0, L > 0 # and chaotic behavior is indicated. 
Downwelling draws in magnetic field-line footpoints; basins of 
downwell ing are field-line footpoint attractors. For upwelling 3^.u < 0, 

L < 0, and 'normal' behavior is indicated. Upwelling blows away field 
line footpoints; domes of upwelling are field-line footpoint repulsors 
(VOORHIES , 1987e) . L is zero for purely toroidal flow. 

Downwelling at a stagnation point (v g = 0) causes exponential 
growth of B r with time constant [9 r u] - ^; a strongly magnetized unipolar 
region called a core spot (Benton & Voorhies, 1981 discussion) or a flux 
bundle or spot (BLOXHAM & GUBB I NS , 1985) may form in an area of 
downwelling. Upwelling at a stagnation point causes exponential decay 
of B ; a weakly magnetized region may form in an area of upwelling. 
Elsewhere the evolution of B r is complicated by lateral advection: 
laterally variable downwelling may strengthen horizontal gradients V s B r 
whose advection by v s can lead to chaotic field behavior. 

If v is a rigid rotation (e.g., bulk westward flow), then B r (b,t) 
is periodic in time. Between extremes is the quas i -per iodic regime 
characterized by weak meridional shear, non-zonal flow components, and 
regions of up- and downwelling embedded in a bulk flow. For example, an 
existing core spot advected through a mild upwelling spreads out and 
weakens, then intensifies upon encountering a mild downwelling and, if 
not absorbed or disrupted downstream, eventually circumnavigates the CMB 
to reencounter the same regions of up- and downwelling. Finally, 
advection of very steep lateral gradients in Bj- produced by vigorous 
downwelling plumes can yield chaotic field behavior. The steady 
surficial motional induction problem is apparently a fine, easily 
visualized example of deterministic chaos which, depending on , 
exhibits : extreme sensitivity to initial magnetic conditions; periodic, 
qua s i - p e r i od i c , and/or non-periodic behavior; fixed attractive or 
repulsive stagnation points with associated basins of attraction or 
repulsion; and regionally chaotic behavior (VOORHIES, 1988b). 

4 . SOME EFFECTS OF PERSISTENT, SURF IC I ALLY GEOSTROPHIC FLOW 

An interesting conceptual view of secular change results from 
supposing core surface flow is both steady and surficially geostrophic. 
With 3 r u = -vtanO/b, chaos may rule in regions of poleward flow. 
Because downwelling, which may form core spots, is accompanied by 
poleward flow, I expect poleward drifting core spots. Conversely, 
regions of low radial flux density formed by upwelling shift towards 
lower latitudes. The implied flux partitioning mechanism can create an 
axial dipole moment from a non-dipole field with a non -axisymme tr ic 
component; moreover, persistent, surficially geostrophic motional 
induction can cause this axial dipole to grow and fluctuate. 

To see how persistent, surficially geostrophic flow can create and 
fortify an axial dipole moment, note that the time rate of change of the 
axial dipole coefficient for a conventional Schmidt-normalized spherical 
harmonic expansion of the scalar geomagnetic potential near Earth s 
surface (r = a = 6.3712 Mm) is 

3 t g^ = ^<9 t B r (a ( t)cos0>. (20a) 
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With either (6) or (13b) just beneath a SFM, 



3b 


2a~ 


-< a t B r 


(b, t ) cos9> 


-3b 


3 


— <V s *(B r v s cos0) 
2a 3 


+ B r vsin9/b> 


-3b 


3 


2a 


-<B r .vsin0/b>. 
3 r 


(20b) 


Frozen-flux variations in the axial dipole require net transport of 
radial field towards the rotation axis, hence poleward flow. 
Surf icially geostrophic poleward flow requires downwelling (12), so 

0 3b 3 

o t g^ = <B r cos0d r u>. (20c) 

2a 3 

The time rate of change of the axial dipole is proportional to the mean 
product of the downwelling and the axial projection of the radial field. 

Pressure perturbations are single valued, so the horizontal 
components of (11) imply that d r u has no axisymmetric component; 
therefore, only the non - axi symme t r ic part of the radial field B r na 
contributes to (20c) for tangentially geostrophic flow. Tangentially 
geostrophic downwelling where B r na cos0 ^ 0 (or > 0) decreases (or 
increases) This can create a normal (or reversed) polarity axial 
dipole from a purely non-dipolar field with a non -axi symme tr i c 
component. This axial dipole will grow until the correlation between 
B r na cos0 and the downwelling changes sign. For steady flow this could 
be quite some time-roughly half the time required for B r na to change 
sign. Indeed, then 5 t 2 g^ = 3b^<B r v s «V s 0 r ucos0>/2a^ . 

It is suggested that surficially geostrophic downwelling in a like 
polarity core spot intensifies the spot and shifts it poleward, 
strengthening the axial dipole; upwelling in a like polarity core spot 
disperses its flux and shifts it equatorward, weakening the axial 
dipole. Conversely, downwelling in a reversed flux patch intensifies it 
and shifts it poleward, weakening the axial dipole; upwelling in a 
reversed flux patch disperses it and shifts it equatorward, 
strengthening the axial dipole (VOORHIES, 1987d) . If such spots or 
patches drift westward, alternately encountering regions of steady (or 
merely persistent) upwelling and downwelling, then quasi-periodic dipole 
oscillations are expected. 

Because there is no downwelling at the poles, a core spot should 
neither form at nor reach the highest lat itudes-which should then be 
regions of weak radial field. Because downwelling implies poleward 
flow, there can be no downwelling at a stagnation point-except at the 
equator. Runaway intensification of a core spot is therefore limited by 
the tendency for the flux to shift poleward (towards regions of weaker 
downwelling). Downwelling at an equatorial stagnation point appears 
problematic: though flux advection from adjacent latitudes vanishes with 
v, an equatorial 'dot' might develop from weak equatorial flux density 
after a few times [3 r u]~ . Yet such a point feature can develop no 
absolute flux linkage as it has no area, its source region (the equator) 
has no area, and downwelling just off the equator is accompanied by the 
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usual poleward flow. Moreover, Lorentz forces and magnetic diffusion 
likely become important near such a dot in violation of (12) and (3). 

Suppose a basically non -axi symmetric configuration of radial field, 
perhaps established by a previous episode of convective/diffusive flux 
expulsion, is subjected to a persistent, surficially geostrophic flow 
which is partly non-zonal. The expected value of (20c) may be zero, but 
the chance of (20c) being exactly zero (of there being no correlation 
between axial field B r cos0 of one polarity and the downwelling) is nil. 
So suppose areas of downwelling correlate initially with areas of normal 
polarity. Then normal polarity flux is attracted by basins of 
downwelling, tending to form core spots, but is also shifted poleward, 
forming high-latitude core spots. The poleward shift inhibits runaway 
intensification of core spots and thus runaway growth of the mean square 
radial field. It seems unlikely that all the flux would end up in an 
axisymmetric configuration; indeed, as normal polarity field-line 
footpoints become stranded poleward of the downwelling foci, dipole 
growth should slow. It seems more likely that flux partitioning would 
continue until some maximum eligible fraction of the normal polarity 
flux has been shifted poleward. Roughly equal partitioning of the mean 
square radial field between dipole and non-dipole configurations might 
be a preferred, possibly metastable state; however, this has not yet 
been demonstrated and continued downwelling might transport reversed 
flux poleward. Perhaps more importantly, nearly zonal advection of 
high- lat i tude core spots through the high- lat i tude fringes of 
alternating regions of up- or downwelling would induce quasi-periodic 
dipole oscillations. With a few core spots and a few foci of up- and 
downwelling, local effects may dominate dipole oscillations at the CMB. 
Yet the global dipole oscillations would appear relatively more 
prominent at Earth's surface. 

This simple picture can be enriched by magnetic flux expulsion: the 
entrainment of toroidally magnetized fluid into the magnetic boundary 
layer, toroidal to poloidal field conversion by the implied laterally 
heterogeneous vertical motion, and diffusion across the CMB. Weak field 
flux expulsion might occur at the low latitudes favored by surficially 
geostrophic upwelling (estimates mapped and posted by VOORHIES (1988a) 
show foci of steady surficially geostrophic up- and downwelling to be 
confined to within about +.3 0° of the equator) and perhaps at the high 
latitudes seemingly favored by density and angular momentum poor plumes 
from the inner core boundary. So the simple scenario should be modified 
to include radial flux diffusion d r (b,t) (which may also help regulate 
core spot amplitudes by eliminating strong curvature r V rB r ) and 
effects of the strong field with appreciable curl expected at depth. 
Sufficiently vigorous vertical advection of such fields may generate 
Lorentz forces that violate even surf icial geostrophy and may thus 
change such a flow. Yet such changes may but punctuate intervals of 
ordinary convective vigor when the Lorentz forces implied by warping a 
strong toroidal field tend to oppose the driving vertical motion. Then 
the simple scenario may commonly hold to a fair approximation. 

To help complete the simple picture, suppose that despite such 
opposing forces, vertical motion at depth grows locally strong enough to 
not only thin the magnetic boundary layer, but to bring toroidally 
magnetized fluid close to the CMB at some intermediate latitude; then 
toroidal to poloidal field conversion by vertical motion and flux 
diffusion can form core spot pairs. If< upon restoration of a 
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surf icial ly geostrophic balance, the spot foci remain in the region of 
upwel 1 ing , then such pairs should weaken and drift equatorward. Though 
such spot pairs are distinguished from the unipolar core spots discussed 
above, the distinction may be artificial. Low-latitude flux expulsion 
and ageostrophic effects may form a spot pair that splits into different 
hemispheres; such spots should be repelled by the upwel ling. Subsequent 
restoration of a surficially geostrophic balance and interactions with 
regions of surficially geostrophic downwelling may lead to further 
separation of the originally paired spots and the formation of seemingly 
unrelated unipolar regions. 

In the long view, a statistically steady, surficially geostrophic 
flow regime characterized by quasi-periodic dipole oscillations could be 
punctuated by vigorous convection (overstability being preferred in 
fluids with small Prandtl and magnetic Prandtl numbers (CHANDRASEKHAR, 
1961) and emergence of one or more bipolar regions. If these spot pairs 
emerge at low latitude and drift into hemispheres of (either like or) 
opposing polarity, subsequent flux partitioning upon restoration of a 
new, persistent, surficially geostrophic flow may dramatically alter the 
dipole moment. Indeed, if enough flux of opposing polarity were 
involved, this could appear as an excursion, or even a reversal, of the 
geomagnetic axial dipole. 

The simple scenario of approximately steady, frozen-flux, 
surficially geostrophic core flow predicts: (1) low-latitude regions of 
low flux (and occasional spot pairs); (2) high- lat i tude unipolar core 
spots; (3) polar regions devoid of spots; and (4) a strong axial dipole. 
Many maps of the broad-scale radial field at the CMB (constructed by 
many workers since BOOKER [1969]) reveal such features. It follows that 
this scenario deserves closer scrutiny. To this end, and in search of a 
quantitatively acceptable explication of recent SV, subsequent papers 
document an effort to describe broad-scale models of the slowly varying 
observed geomagnetic field in terms of (piecewise, statistically) 
steady, optionally surficially geostrophic, core surface motions. 

5. SUMMARY 

The simple source-free mantle/ frozen-flux core (SFM/FFC) earth 
model has been used to review and extend the kinematic and elementary 
dynamic theory of geomagnetic secular change. The ROBERTS & SCOTT 
[1965] equation was used to prove that the mean square radial magnetic 
flux density averaged over a FFC , <B r (b,t) 2 >, changes if and only if 
fluid downwelling, -V s «v g (b,t) = 3 r u(b,t), correlates with squared 
radial magnetic flux density, and is thus conserved by any purely 
toroidal FFC surface flow. For a SFM, the contribution from the broad- 
scale portion of the radial magnetic flux density to <B r (b,t) 2 > can be 
easily estimated from broad-scale models of the evolving geomagnetic 
field and used to check the compatibility of such models with the no- 
upwelling hypothesis. 

The kinematical forward problem posed by steady motional induction 
at the top of a FFC was solved analytically (albeit not in closed form) . 
The sign of the Liapunov exponent describing evolution of the radial 
magnetic flux density B r (b,t) was shown to be the sign of the 
downwelling d r u (b) -raising the possibility of steady motional induction 
of geomagnetic chaos. Steady motional induction at the top of a FFC was 
found to a be fine, easily visualized example of deterministic chaos 
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which, depending on v g (b), exhibits extreme sensitivity to initial 
geomagnetic conditions; periodic, quas i -periodic , and non-periodic field 
behavior; fixed basins of magnetic field-line footpoint attraction or 
repulsion surrounding fixed stagnation points; and regionally chaotic 
field behavior. 

Downwelling implies poleward flow for surficially geostrophic core 
flow. Some aspects of persistent, surficially geostrophic core motions 
were described, including the flux partitioning mechanism whereby the 
formation and poleward drift of core spots in regions of fluid 
downwelling can create and fortify an axial dipole from a field that is 
at least partly, and could be entirely, non-axisymmetr ic . Features 
expected to result from persistent, surficially geostrophic flow are 
evident in many mapped estimates of the broad-scale radial field near 
the top of the core, and may be present in records of pa 1 eomagne t ic 
secular variation. However, magnetic flux diffusion, Lorentz forces, 
and t ime -dependent flow are expected to play key roles in the correct 
interpretation of very long term SV and geomagnetic reversals. 

In one long-term scenario, intervals of persistent, partly non- 
axisymmet r ic , surficially geostrophic core flow characterized by quasi- 
periodic global dipole oscillations and regionally chaotic field 
behavior are intermittently punctuated by major changes of the flow in 
the outer core. The instability envisioned is of internal origin-a 
natural consequence of thermo-compositional convection within a rapidly 
rotating, roughly spherical annulus of strongly magnetized liquid metal- 
alloy subject to very slowly changing, inhomogeneous cooling from above 
(due to thermal instability of the deep mantle) and a very slowly 
changing, heterogeneous flux of buoyant material from below (due to 
condensation of the inner core) . It is characterized by uncommonly 
strong poloidal flow within the upper core, entrainment of toroidally 
magnetized fluid into the magnetic boundary layer (say the upper 90 km 
of the core), conversion of toroidal to poloidal field by the vertical 
motion, flux diffusion across the CMB , and formation of core spot pairs 
-particularly, but by no means entirely, at the low latitudes favored by 
surficially geostrophic flow. There may be many such pairs which split 
ageos trophica 1 ly ; if enough spots of sufficient flux happen to drift 
towards, if not into, hemispheres of opposing polarity, then the process 
appears to be either an excursion, or a reversal of the axial dipole. 
Relaxation to a new pattern of persistent, partly non-axisymmetr ic , 
surficially geostrophic core flow leads to intensification and poleward 
drift of such spots, formation of unipolar core spots at high latitudes 
and reestablishment of an axial dipole. The polarity of the new dipole 
is the polarity of the spots which happen to correlate with the new, 
surficially geostrophic pattern of downwelling. 

This is of course but one of many conceivable scenarios in need of 
development and, more importantly, testing. Many paleomagnet ic data and 
geomagnetic data can be brought to bear upon speculations, properly 
framed hypotheses about the core, and geodynamo theory in general if the 
physical suppositions, mathematical means, and numerical methods needed 
to connect hypotheses and observations are developed. In order to test 
simple hypotheses such as (statistically) steady, piecewise steady, or 
persistent surficially geostrophic flow and, perhaps more importantly, 
in search of a quantitatively acceptable explication of recent SV, a 
method has been developed to connect the SFM/FFC kinematics of secular 
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change to broad-scale models of the evolving geomagnetic field. This 
method is described in the next paper. 
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APPENDIX A 


The magnetic induction equation, though familiar to some, is not as 
widely known as can be hoped. Moreover, some textbook derivations of it 
(e.g., ROBERTS, 1967; JACKSON, 1975; GUBBINS & ROBERTS, 1987) skip a few 
points of interest in geomagnetism; in particular, a common rationale 
for using Ampere's law is inadequate for some geomagnetic purposes. So 
we return to the differential form of the macroscopic Maxwell equations. 

Let q represent the electric charge density scalar (with dimensions 
of C/m 3 ), E the electric field vector (V/m) , D the electric displacement 
vector (C/m 2 ), J the electric current density vector (A/m 2 ), B the 
magnetic flux density pseudo-vector (T) , and H the magnetic field 
strength pseudo-vector (A/m) . The Maxwell equations describing the 
position r and time t -dependent electromagnetic field are then 

V#D = q (Ala) V®B = 0 (Alb) 

VxE = -3 t B (Ale) VxH = J + 3 t D (Aid) 

where 3 denotes partial differentiation with respect to the ensuing 
subscripted variable, V* the divergence operator, and Vx the curl 
operator. The divergence of the Ampere -Maxwe 1 1 law (Aid) is charge 
conservat ion 

9 t q + V*J = 0. (Ale) 

In the reference frame K' moving with a physical medium the 
macroscopic constitutive relations are 

B' = | 1 *H' (A2a) D' = E*E' (A2b) J ' = a® E' (A2c) 

where y(r',t') represents the second-rank magnetic permeability tensor 
(H/m) , (r ' , t ' ) the dielectric permittivity tensor (F/m) , and g(r',t') 
the electrical conductivity tensor (S/m) of the medium. In (A2) the dot 
product • indicates the first-rank result of a second-rank tensor 
operating on a first rank vector (or pseudo-vector). For non-linear 
media \i, and a depend implicitly on H' and E' . For steady media y, 
and c are independent of t'. For homogeneous media \X, and Q are 
Independent of r' . Because H ' and E' are generally functions of both r' 
and t' a medium which is steady, or homogeneous, or both (i.e., uniform) 
will generally be linear. For isotropic media y, £., and Q are diagonal 
matricies, each with three identical elements, and are described by the 
scalars u, e, and o. In an isotropic medium, the speed of light is c ^ 
(|l£) , the charge relaxation time is s e/a, and the magnetic 

diffusivity is T| = (|ia)” 3 . 

Jump conditions across the interface between two different media 
are derived from the integral forms of (Al) (see, e.g., JACKSON, 1975). 
Let n be the unit vector normal to the interface, A n s A»n, and A s s A - 
A n n. Then these conditions are continuity of B n and E s and jumps {D n ) = 
X and { H g } = Cxn, where X is the idealized surface charge at the 

interface and C is the idealized surface current at the interface. 
Given \± , 0, and time-dependent boundary conditions, suitable initial 

conditions at t' =0 for solving (Al, A2 ) are E'(r f ,0) and B'(r',0). 


17 



For a medium moving with uniform velocity u relative to an inertial 
reference frame, the comoving frame K ' is inertial and equations (Al) 
and (A2) are implicitly covariant under the special relativistic Lorentz 
transformation (see, e.g., JACKSON, 1975). For the outer core, the 
fluid velocity v is not uniform so K' is not inertial. Moreover, v is 
defined relative to the non-inertial reference frame K fixed to the 
imperfectly rigid, solid Earth which is in non-uniform rotational and 
orbital motion relative to an inertial reference frame. 

The portions of the electromagnetic field of geophysical interest 
are described qualitatively by their length scale of spatial variation 
X, time scale of temporal variation X, and typical speed U = X/X ~ |v| . 
Slow changes of the geomagnetic field attributed to core sources are 
described by X < 2nb - 22 Mm and x > 1 year = 32 Ms, so U < 0.69 m/s. 
Typically X > X^ « 2rcb/14 = 1.6 Mm, T > X^ ~ 10^ years, and U ~ 5x10’^ 

u u 4 

m/s. Clearly U is much less than the speed of sound c 0 ~ 10^* m/s within 

o * 

the Earth and the vacuum speed of light c Q = 3.0x10 m/s; however, Earth 
is not vacuum, so geophysical scale analyses comparing U with c Q may be 
irrelevant. If e o < e < 10“^ F/m and if fi = (1 0 = 4 tcx 10~^ H/m, then c > 
9xl0 3 m/s >> U; but do we know e? Because hydrodynamic stress and key 
bulk material properties (elastic moduli) used to describe sound are 
macroscopic representations of microscopic electromagnetic interactions 
between particles, sound is an electromagnetic phenomenon. So the top 
speed at which information can propagate electromagnet ical ly within the 
Earth is not less than the sound speed. Then c > c s >> U and progress 
can be made even in the absence of measured ultra-low frequency values 
for E and |l of deep-Earth materials. 

With U << c < c Q the general transformation from K' to K should 
approach the Lorentz transformation which, in turn, approaches the 
Galilean transformation B' = B and E' = E + vxB; moreover, J' = J - qv 
and charge invariance imply q' = q. Scale analysis of Faraday's law 
(Ale) yields |e| ~ X\b\/X = u| b| , so the Ampere -Maxwe 1 1 law (Aid) is 
approximated by Ampere's law (\^<H = J) because the displacement current 
is allegedly small compared with the curl of the field strength: 
la t D| /I VxhI - e|BlnX/|B|x - (U/c) 2 « 1 (ROBERTS, 1967 ). This should 
hold even with c << c Q ; yet |3 ^d|/|Vxh| = 1 in regions where J = 0 and 

the geomagnetic field is supposedly scaloidal . Conditions under which 
Ampere's law is useful are offered in Appendix B: for steady e and a it 
is shown that x must be very much greater than Tg for this ultra-low 
frequency, electromagnet ically quasi-steady approximation to be useful. 

For electrically neutral media q = 0 and J' = J - qv = J. More 
generally, lqv| = |(V*eE)v| - eIeIu/X, - U^IbIe/A, and scale analysis of 
Ampere's law yields | b| / jxA - |j|; therefore, |qv|/|j| - U^/c^ and J' is 
approximated by J (ROBERTS, 1967). This holds even with c << c Q . 

These approximations yield the pre-Maxwell equations appropriate to 
the magnetohydrodynamics of an isotropic medium 

V*EE = q (A3 a) V*B = 0 (A3b) 

VxE = -a t B (A3c) 
and Ohm's law 

J = a (E + vxB) . ( A3e) 


Vx (B/(l) = J (A3d) 
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The vxB term is an order-zero relativistic effect much like classical 
kinetic energy m 0 v 2 /2. (Recall that for invariant rest mass m Q , E = 
m 0 c 2 (1 -v2 /c 2 ) -l/2 0 = mQC 2 + / 2 + terms of order v 2 /c 2 (OKUN, 1989)}. 

Jump conditions across the fixed interface between two different 
media in which (A3a-d) hold are derived from (A3b), (A3c), the 

divergence of (A3d) , and (A3d)^. These conditions are continuity of B n , 
E s' J n # and ^ uinp = Cxn across the interface. If J is finite, 

then {B /(!} =0. If a is finite at the interface, then finite E, v, and 
B in ( A3 e ) imply finite J; then it seems C = 0 and B s /[1 is continuous. 
If the finitude of G is omitted (as in frozen-flux) , then C ^ 0 and 
discontinuous B g /jl are allowed; then B g may jump even if |1 does not. 
The appropriate initial condition at time t = 0 is B(r,0) . The initial 
condition on E is not important in the quasi-steady approximation and 
must be sacrificed because (A3d) is a singular perturbation of (Aid) . 
Given solenoidal B (A3b) , Ampere's law (A3d) determines solenoidal J; 
further given v, Ohm's law (A3e) determines E. Then Gauss's law (A3a) 
determines q and (eE n ) = I. Comparison of V»J = 0 with the true charge 
conservation law (Ale) confirms that sources of displacement current - 
3j.q (along with qv, and Xv ) are omitted in the quasi-steady 

approximat ion . 

Evaluation of 3 t q from (A3a-e) yields 
a t q = a t V*(eE) = a t V.(x q J - evxB) = 3 t (Vx q »J) - V*3 t (evxB) (A3f) 

= V»3 t (x q VxH - c" 2 vxh) 

which must be negligible compared with q/X q , |j|/X., and Isl/pX. 2 for 
Ampere's law to be useful in conjunction with (A3e) . If the quasi- 
steady approximation is to be strictly compatible with true charge 
conservation (Ale) # then (A3f) must be identically zero. For uniform 
media this would require eV*3 t (vxB) = e9 t (B®o - ^lv»J) = 0, where CO = Vxv 
is the vorticity. This requirement is always met in the comoving frame 
(v = 0); in the limit as e, and thus and c 2 , approach zero; in 
electrically neutral media; and when (B*C0 - is steady (e.g. r 

irrotational sub-relativistic motion of an electrical insulator). The 
unimportance of (A3f) is necessary, but not sufficient, for use of the 
quasi -steady approximation. 

Elimination of E from (A3c) using (A3e) yields 

3 t B = Vx(vxb) - Vx(j/o) . (A4) 

Substitution of Ampere's law (A3d) into (A4) yields the magnetic 
induction equation (la). 

Omission of qv does not imply J = 0. Charge density q is a 
macroscopic average over the microscopic charge carriers of their 
individual charge alone; charge current density J is an average over the 
charge carriers of their individual charge multiplied by their 
individual velocity. Similarly, mass density p is a macroscopic average 
over the microscopic mass carriers of their individual mass alone; mass 
current density pv is an average over the mass carriers of their 
individual mass multiplied by their individual velocity. Because pv is 
dominated by ponderous mass carriers (neutral molecules, 'holes', and 
ions) while J is dominated by mobile charge carriers (conduction and 
free electrons), pv need not be obviously related to J. The electro- 
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magnetic quasi-steady approximation ( V« J = 0) does not imply the 
acoustic quasi-steady approximation (V«pv « V®p Q v = where p o (r) is 
the mean density stratification p o (r) » p(r,t) - p Q (r) = p'(r,t)); 
moreover, a steady mass current density (3 t [pv] = fP 0 v ^ = Po^t v = 
does not imply a steady charge current density (d t J = 0) . 

APPENDIX B 

For media in which £ and a are non-trivial, isotropic, and uniform 
equations (Ala), (Ale), (A2b) , and (A2c) in comoving frame K' require 

eV'»3 t ,E' + oV'*E' = 3 t , q ' + X q _1 q' = 0 (Bl) 

so 

-t'/Tq -t'/T 

q'(r',t') =q'(r',0)e = eV ' «E 7 (r ' , t ' ) = eV'*E' (r ' 0 ) e . (B2) 

Provided X q > 0, both q' and V'*E' decay exponentially towards zero with 
time constant x q = e/G. (Note e, (1, and a are regarded as purely real 
for simplicity) . 

If e and a are non-trivial, isotropic, and steady, but possibly 
inhomogeneous, then in K' , equations (Aid), (A2b) , and (A2c) require 

9 t ,E r + T q ~ 1 E ' = £ _1 V'XH'. ( B3 ) 

Because T q (r')= e(r')/o(r') is also isotropic and steady, the solution 
of this first-order equation is 


E' (r',t') 


-t '/T r 


E' (r' , 0)e 


+ 


t' "(t' - X)/X 

e(r') 1 J V ' xH ' ( r ' , X ) [ e ]dx. (B4) 

0 

Suppose x is positive def inite-which excludes ' ant i-f erroelectrics ' or 
'anti-conauctors . ' Then E' forgets its initial condition at t'=0 after 
a few x q and is effectively the exponentially weighted running average 
of e -1 V'xH' thereafter (BACKUS, 1982). Time integration of (Ale) after 
using (B4) to eliminate E'(r',t'>0) yields 

-t'/T 

B' (r ' , t ' ) = B(r' , 0) - V'xE' (r' , 0)T [1 - e ] - 

t' T -(T - T)/T 

V ' xe { r ' ) " 1 | {1 V ' xH ' ( r ' , X ) [ e ]dx}dT. (B5) 

0 0 

If part of H' does not vary exponentially on time scales near t q 
when t' >> X q > 0, then only values of this part of e -1 V'XH' within a 
few X q of t' contribute much to E' because the weight factor in (B4) 
definitely decays exponentially. Slow or secular variations in E' and 
H' on time scales long compared with x q may then be considered quasi- 
steady: E ' « X q e -1 V'XH / or, by (A2c) , 
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VXH' (r ' , t ' ) = J' (r',t') ^ DD ' 

which is Ampere's law. Physically, the Maxwell displacement current in 
(Aid) is omitted in (B6) because fast changes in D' , hence in q', E', 
and thus B' by (B5), contribute negligibly to the slowly varying 
portions of interest. A more rigorous derivation of (B6) follows. 

We are concerned with the comoving frame K' but omit the prime 
notation for now. Again suppose e and O are non-trivial, isotropic, and 
steady; suppose H(r,t) is also real. During the time interval from 0 to 
t, H(r,T) is the sum of a constant H 0 (r), a linear trend ii 0 (r) [x-t Q ] , 
and a fluctuating portion h(r,T) = 5^ k k (r) sin (0) k I+<|> k ) : 


H(r,T) = H 0 (r) + H c (r) [T 


t c ] + £ h k (r) sin(0) k x + <|> k ) . 

k=l 


(B7) 


If (Oi, were 2Jlk/t, then the sum would be the Fourier representation of 
the zero mean detrended magnetic field strength during the interval . We 
do not insist on this choice so as to allow flexibility in defining the 
mean, the trend, and any fundamental period; indeed, a transform over 
continuous CO may replace the discrete series. With (B7) the three 
contributions to the integral on the right of (B4) are: 


t - ( t-x) /x, 

I VxH 0 (r)e 
0 


-t/Xg 

q dx = VxH 0 (r)X q [l-e ] 


t . 


- ( t-X) /X, 


J VxH (r) [x-t ] e 


-t/tq 

q dx = VxH 0 (r)X q {t - [t 0 +T q ][l-e ]) 


( B8a) 


(B8b) 


t - ( t -X ) /X q 

j Vxh k (r) sin((0 k x + <|> k )e 
0 


= Vxh k (r)X q [l+(0 k 2 Xq 2 ) 1 (sin (CO k t+<|> k ) - ® k X q cos (0> k t+<t> k ) 

_t/ ^q 

- [ sin<)> k - (O k XqCos<t> k ]e ) (B8c) 


ty ^q 

= Vxh k (r)X q (l+(O k 2 Xq 2 ]" 1 (sin(t) k - (O k x q cos<t) k ] [1 - e ] 

where the last step would follow from (B8c) only if the 0^ were 2rck/t . 
With (B8a-c) , multiplication of (B4) by a(r) yields 


J(r, t) 


-t/x 


q 


J(r, 0)e 
+ VxH 0 (r) {t 


-t/X, 


+ VxH 0 (r) [1-e 


-t/x. 


[t 0 +Xq] [ 1 e 


] 

]) 
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CO 


+ Z Vxh k (r) [l+ 0 ) k 2 t q 2 ] 1 {sin(co k t+<|> k ) - w k x q (cosco k t+<|> k ) 
k = 1 

-t/T q 

- [sin(j) k - 0) k T q cos<|) k ]e }. (B9) 

When t » X q , exp(-t/X q ) » 0 and 

J(r,t) - VxH 0 (r) + VxH Q (r)[t - t Q - x q ] + 

00 

1 Vxh k (r) [l+(0 k 2 x q 2 ] _1 [sin(C0 k t+()) k ) - w k x cos (co k t+<|> k ) ] . (BlOa) 
k=l 

For t >> X q , one may choose t Q : X q << It - t Q | . So choosing implies 

~ “ [t - t Q ] . If Ico^liq << 1 for all non-trivial Vxh k , then 

to order zero in the small quantity co k x q [ 1+C0 k 2 x q 2 ] _1 « co k X q « 0, 

. CO 

J(r,t) * VxH Q (r) + VxH 0 (r) [t - t Q ] + X Vxh k (r) sin (CD k t+$ k ) 

k=l 

= VxH(r, t ) . (BlOb) 

Restoring the primes converts (BlOb) into Ampere's law in the comoving 
frame (B6) . Ampere's law holds only for frequencies |(D k l small compared 
with X q so it is strictly an ultra-low frequency approximation. 

If CD k were 27tk/t, then l© k lx q << 1 would imply x << 27lt/k. 
Ampere's law is thus broken by any non-trivial Vxh k with k > 2rct/x 
If, however, the Vxh k (r) are negligible for k > K, then the sum can §e 
truncated at k = K; then Ampere's law will be useful provided co^t = 
Kx q /2rct << 1 (i.e., t >> Kx q /2 k ) . More generally, to^ {or continuous CO) 
can be very high; yet Ampere's law will hold if the ^<h k are negligible 
for I C0 k I » Tq' 1 , Clearly the use of Ampere's law is tantamount to 
truncating the sum over frequencies; higher frequencies require (Aid). 

When the total electromagnetic field has finite energy in finite 
volume , it has finite energy density at all but a finite number of 
singular point sources with finite charge and finite magnetic moment. 
Then we may insist that the macroscopic averaging procedure give h k 2 (jr) 
which approach zero as IcOjJ approaches infinity. Unfortunately, it is 
not clear that the cutoff frequency above which the h k are negligible 
will be small compared with X q . On Earth's surface, high-frequency 
electromagnetic oscillations comprising the solar and geothermal fluxes 
break Ampere's law, but have far less energy density than the main 
geomagnetic field; yet within Earth's core, the electromagnetic energy 
density associated with high-f requency inter-molecular collisions (hence 
hydrostatic pressure) vastly exceeds that of the main field. 
Nevertheless, when £ and 0 (and )i) are steady, the temporal linearity of 
the Maxwell equations in the comoving frame ensures that different 
frequencies are linearly independent. Then it is useful geomagnet ically 
to filter out high-frequency electromagnetic oscillations by truncating 
the sum in (BlOb). This is accomplished experimentally either by using 
low-pass magnetometers or by averaging the results in time. 
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Because Ampere's law holds to order zero in the relaxation time Tq, 
it is tempting to think of it as the limit as E and thus c approach 
zero; then (A3f) would approach zero. Although the quasi-steady 
approximation (A3d) and (A3e) follows from treating c and E (but not 
|l) as if they were zero, it does not require E to be treated as if it 
were zero in (Ala) or (A2b) . Use of Ampere's law merely requires 
| a t D| / I j| = I 3 t (eE) I / I CeI ~ I xv»3 t (£B) I / I OBl ~ T q /T << 1 ■ Galilean 
invariance, Ohm's law, and the omission of qv, hence the quasi-steady 
approximation, further require terms of order U 2 /c 2 to be negligible. 

In the special case when o is identically zero in some region, then 
J is zero, q is constant, and T q is infinite. Then t cannot be greater 
than T - and Ampere's law holds only if = 0. Yet Ampere's law can 
still be used if displacement current sources outside the region, which 
give rise to solenoidal 3^_D = Vxh within the region, contribute 
negligibly to the portion of H of interest. Then this portion of H is 
irrotational and originates outside the region of interest. 

Curiously, one rationale for the omission of displacement currents 
( I a t D| / 1 VxhI - (X/T C) 2 = U 2 /c 2 « 1) typically fails in regions where H 
is irrotational (where \i is homogeneous and B is scaloidal). Moreover, 
this rationale is often used to justify the use of Ampere's law in 
conductors, where it is not sufficient. Fortunately, conditions under 
which U 2 /c 2 << 1 but t * T q appear to be but rarely encountered in solid 
Earth geomagnetism; however, they appear to be common where fluctuations 
in the conductivity and dielectric properties of sub- r e la t ivi s t ic 
plasmas are of considerable interest. Such non-linear media are 
considered briefly in Appendix C— as are steady, anisotropic media. 

APPENDIX C 

Real media are at best statistically steady and statistically 
isotropic, so it seems worth considering the case of inhomogeneous, 
anisotropic, and unsteady media. This includes non-linear media in 
which the dependence of , G, and y on electromagnetic fields is 
implicit in the dependence of their elements on position r and time t. 
The Ampere -Maxwe 1 1 law (Aid) and the anisotropic, possibly non-linear, 
constitutive relations in the comoving frame require 

e*5 t E + (a + a t e)*E = Vxh = Vx[^ _1 *B] (Cl) 

where the prime notation is again supressed and only media in which 
possibly non-linear second-rank operators E. ^ r G ^ and y exist are 
considered. Define the tensor relaxation time operator t q (r,t) via 

J[e -1 *(a + 3 t e)]dt = tx q _1 (C2) 

where the integrand is the second-rank result of matrix multiplication. 
With (C2 ) and (A2c) , the solution of the first-order equation (Cl) is 
written 

-tx - 1 
tl q 

E (r, t ) = {e }»E (r, 0 ) + 
t -(t - T)T _1 

J {e “ )«{e(r,X) _1 «Vx[(i.(r,X) i *B(r,X) ] }dx (C3) 
0 = 
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where expi-ti^' 1 ) is itself a matrix operator that varies with (r,t). 

If it is possible to represent the elements of Iq” 1 / , y" 1 , and 
B in terms of their mean values, trends, and fluctuating portions during 
the interval from 0 to t, then one may proceed by analogy with Appendix 
B. Substitution of these expansions into (C3) reveals that the curly 
bracketed term on the far right of (C3) 


U(r,T) 2 {e (r, t) 1 *Vx[n(r,x) 1 »B(r,T)]) 


(C4 ) 


will include terms proportional to: (a) (T-t Q ) n for n = 0, 1, 2, or 3; 

(b) the products of single oscillations and (t-t 0 ) rn with m = 0, 1, or 2; 

(c) the products of two oscillations and (x-t^ 1 with i = 0 or 1; and 

(d) the products of three oscillations. The interactions of £., y, and B 

generally will contribute mean, secular, and oscillatory terms to the 
expansion for U. For example, under resonant conditions in which the 
two oscillations in the i = 0 terms of case (c) have the same frequency 
and phase, there generally will be a non-zero contribution to the mean 
value of U. Therefore the mean value of U will generally not be the 
same as obtained by combining the mean values of t, y, and B according 
to (C4) . This can be the case even if is steady, provided y and B 

are not. And the exponential operator in (C3) must contain the trend 
and oscillations in as well as its mean value. 

It therefore seems extremely unlikely that Ampere's law will hold 
for unsteady media; however, the problem does not arise for steady media 
even if they are anisotropic. And from the macroscopic perspective, 
unsteady media often appear to be near a statistically steady state. If 
this is the case, then the macroscopic constitutive relations (A2) can 
be considered definitions of the statistically steady, macroscopic 
properties of the medium. 

Of particular interest in solid-Earth geomagnetism are cases in 
which £ and g in (A2b) and (A2c) are taken to be (statistically) steady, 
albeit possible anisotropic and inhomogeneous. To the extent that £ and 
g are steady, 1 ' is steady. If H is used instead of y*B, then 
one can proceed by analogy with Appendix B, albeit in matrix notation, 
without placing restrictions on y. Indeed, with possibly complex £_ (r) 
and g(r) , and with 

i« k T 

H (r , t) = H 0 (r) + H Q (r) [x - t Q ] + I h k (r)e 

k 

equation (C3) becomes 


-tx. 


-1 


-tx. 


-1 


E(r,t) = 


_q w lq -i 

{e " }»E(r,0)+[I-e “ 


-tx. 


-1 


+ { [ t - t Q ]I - l q + [e 


JTtoI - }*C 1VxH o 


-tx, 


lq 1 , 

+ [I - e ] • [T_ 1 + i(D k I]«E i «Vx[£ h w 


ico kt 


_q 


(C5) 


where I is the identity matrix. If, for large t, exp(-tt q 1 ) = Q, then 
(C5) simplifies to 


1 




f 


i 


S 
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(C6 ) 


E(r f t) ~ o 1 #Vxh o + {[t - t 0 ]l - T q }*a 1 Vxh q 

iCD^t 

+ [t + ico^I] Vxh^e 

= = = k 

If the frequencies of interest are so. JLow that [£ q ^ + iCO^I] ^ and 

if the elements of possibly complex are negligible compared with [t 

- t Q ]I, then (C6) reduces to Ampere's law. In such cases, Ampere's law 
is a justifiable ultra-low frequency approximation to the Ampere -Maxwell 
law for inhomogeneous, anisotropic media with unsteady magnetic 
permeabi 1 i ty . 

The possibility of unsteady y yields a flexibility which may be 
useful when dealing with time-dependent crustal rock magnetization 
M (r , t ) , where M(r,t) = [ 1^ 0 ” (r , t ) - I]*H(r,t). Note that = 

Zi^i( r ,t) can be a sum of matricies. One such portion could be the 
product of a matrix @(r,t) and the diagonal matrix with elements 
inversely proportional to H(r,t). This would yield a hard magnetization 
which is independent of the 'inducing' H, but, depending on the form of 
@, might depend on other macroscopic parameters such as temperature 
T(r,t). Other portions of y(r,t) can be constructed to mimic other 
sorts of highly non-linear behavior, including thermo- viscous remanent 
magnetization . 
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